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We investigate two schemes for pumping spin adiabatically from a ferromagnet through an inter- 
acting quantum dot into a normal lead, which exploit the possibility to vary in time the ferromagnet's 
magnetization, either its amplitude or its direction. For this purpose, we extend a diagrammatic 
real-time technique for pumping to situations in which the leads' properties are time dependent. In 
the first scheme, the time-dependent magnetization amplitude is combined with a time-dependent 
level position of the quantum dot to establish both a charge and a spin current. The second scheme 
uses a uniform rotation of the ferromagnet's magnetization direction to generate a pure spin current 
without a charge current. We discuss the influence of an interaction-induced exchange field on the 
pumping characteristics. 

PACS numbers: 73.23.Hk, 72.25.Mk, 85.75.-d 



I. INTRODUCTION 

The last decades have seen intense research in the field 
of spintronics, a branch of physics and electrical engi- 
neering whose aim is to combine electrical and magnetic 
functionalities in the same solid-state system by exploit- 
ing the spin degree of freedom of the charge carriers!^! 
In particular, mechanisms to generate and control spin 
currents have been the object of several theoretical and 
experimental investigations. A subset of these studies, 
which is particularly relevant to the present paper, fo- 
cused on quantum-dot spin valves, that is devices made 
up of a quantum dot sandwiched between two ferromag- 
netic leads. Experimentally, these types of systems have 
been realised with self-assembled In As quantum dotsjpsl 
small metallic grainsptEH nanowires™ nanotubea^HIl] 
and molecular systems.^ Quantum dots contacted with 
ferromagnetic leads have also attracted a considerable 
theoretical interest EsHUl 

A possible mechanism to generate a pure spin current, 
that is a spin current without an associated charge cur- 
rent, relies on a ferromagnet with a rotating magnetiza- 
tion. This mechanism is the corner stone of recent spin- 
battery proposals EStSU in a related study, it has been pre- 
dicted that the rotation of the magnetization direction of 
a magnetic quantum dot that is sandwiched between a 
ferromagnet and a normal lead gives rise to charge pump- 
m g[53] Similarly, spin pumps using time-dependent mag- 
netic fields acting on a molecule have been proposed! 54 * 55 !. 
These time- dependent transport problems can be formu- 
lated within the theo ry of charge and spin pumping in 
mesoscopic structures J 42 l 56 l EH However, pumping usually 
refers to the situation when the properties of the conduc- 
tor and not those of the leads (as it is the case for the 
spin battery) are varied in time. The regime of adiabatic 
pumping is realised when the period of the time variation 
is long compared to the characteristic dwell time of the 



carriers in the system. This regime is particularly inter- 
esting from the theoretical point of view since it leads to 
an understanding of the non-equlibrium caused by the ex- 
plicit time dependence of the system without introducing 
the additional complications associated with higher-order 
terms in the frequency of the time-dependent parameters. 

In this paper we focus on a spin battery realised in 
a structure composed of a quantum dot, tunnel coupled 
to a ferromagnetic and a non-magnetic lead. The spin- 
battery operation is due to the time-dependent magneti- 
zation of the ferromagnetic lead. We employ a diagram- 
matic real-time theory for adiabatic pumping through 
quantum dots with ferromagnetic lead d 42 * 7 ^ that we ex- 
tend to account for time-dependent magnetizations. This 
approach consists in a systematic perturbative expansion 
in powers of the tunnel-coupling strengths and of the 
pumping frequency, while treating the on-site Coulomb 
interaction on the quantum dot exactly. We consider two 
different pumping schemes: First, we choose the ampli- 
tude of the magnetization of the ferromagnetic lead and 
the level position of the dot as pumping parameters. Sec- 
ond, we pump by changing periodically the direction of 
the magnetization. A sketch of the system under inves- 
tigation and of the two pumping schemes is shown in 
Fig.0 

The main results we find are the following. In the 
first pumping scheme the pumped spin current is, in gen- 
eral, accompanied with a finite charge current. There 
are, however, special values of the system parameters, 
for which the pumped charge current turns out to be 
zero due to a cancellation of two counteracting contribu- 
tions to the overall pumped charge. As a consequence, in 
the first pumping scheme the generation of a pure spin 
current requires a fine tuning of the system parameters 
. This is profoundly different to the second pumping 
scheme in which the pumped charge vanishes for any 
choice of the system parameters. We are able to derive 
a compact analytical formula for the pumped spin cur- 
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FIG. 1. (Color online) Sketch of the F-dot-N structure. It 
consists of a quantum dot, tunnel coupled to a ferromagnetic 
and a non-magnetic lead, which are kept at the same chemical 
potential. 

(a) Schematic description of pumping with the amplitude of 
the magnetization of the ferromagnetic lead and the level- 
position e of the dot. 

(b) Schematic description of pumping with the x- and y- 
component of the magnetization of the ferromagnetic lead. 



rent and to discuss its dependence on the gate voltage, 
the ratio of the tunnel couplings to the ferromagnet and 
the normal lead, and the exchange field that acts on the 
quantum dot spin as a consequence of a spin-dependent 
tunnel coupling of the dot level to the ferromagnet. 

This paper is structured as follows. After introducing 
the model of the system in Sec. [TTJ we extend in Sec. [TTT] 
the diagrammatic real-time technique to calculate the 
pumped charge and spin through systems including fer- 
romagnets with time-dependent magnetization. The dia- 



grammatic rules resulting from this are given in Sec. IV 



The main part of the paper is Sec. |V| where the results 
for both pumping scheme are presented. We conclude by 
summarizing our results in Sec. |VI| 

To keep all formulas transparent we set H = 1 through- 
out the paper. 



II. MODEL 

The Hamiltonian of the system reads H = H^ot + 
-ffiead n + -^lcad f + -^tun- The single-level quantum dot 
is described by the Anderson impurity model 



Hdot = ^ e d\ d a + U ivf m , 



(1) 



where e is the spin-degenerate energy level of the dot, 
U the on-site Coulomb interaction, d a {d\) the annihi- 
lation (creation) operator of an electron on the dot with 



spin a =f, 4-, and n a — S a d a the corresponding number 
operator. The Hilbert space for the quantum dot is four- 
dimensional. The basis states are denoted by |0), |t), 
and |d) = d|d||0) corresponding, respectively, to the dot 
being empty, singly occupied with spin a =t, i, or doubly 
occupied. 

The leads are modeled by reservoirs of non-interacting 
electrons (kept at the same chemical potential fi which 
we set to zero) with Hamiltonians 



-fflead N = £fe C' k<y Cka 



Hi 



IcadF 



(2) 
(3) 



k.a 



where Cka ( c htr) 1S ^ nc annihilation (creation) operator 
for an electron with spin a and wave vector k in the 
normal lead, while a^a ( a k a ) is the annihilation (cre- 
ation) operator of an electron with majority /minority 
spin a = ± and wave vector k in the ferromagnetic 
lead. The density of states g^(ui) = J2 k S(cj — e&) 
in the normal lead is independent of spin. The ferro- 
magnet, on the other hand, is described by majority- 
and minority-spin bands, E ka = e*. + E a , that are 
shifted relative to each other by a finite Stoner split- 
ting AE = Ek- — Ek+- This leads to a spin-dependent 
density of states q F -{uS) = Qf,+ (u — AE). The de- 
gree of spin polarization at energy co is characterized by 
P(uj) = [gF )+ (w)-0 P) _(w)]/[0F,+ (w) + £F,-(w)]. We 
remark that the majority and minority spin direction in 
the ferromagnet (denoted by a — ±) may, in general, be 
different from the spin quantization axis (with a =t, 4-) 
that we choose for the quantum dot and the normal lead. 
The corresponding Fermi operators are connected via the 



transformation a| 



with 



e-^/ 2 cos((9/2) e^/ 2 sin(0/2) 
-e- J¥> / 2 sin((9/2) e iv / 2 cos(9/2) 



(4) 



where the polar angle 9 and azimuthal angle ip de- 
fine the ferromagnet 's magnetization direction e p = 
(sin8 cos ip, sin 9 sirup, cos#) T in the (time-independent) 
coordinate system with the z-axis chosen along the spin- 
quantization axis of quantum dot and normal lead. 

The tunnel coupling between dot and leads is described 
by the spin-conserving tunneling Hamiltonian 



fftu„ = (^ N c l<y d ° + Vf a l<y d ° + hj 



(5) 



where Vj$ and Vp are the energy- and spin-independent 
tunnel-matrix elements. Tunneling introduces a finite 
lifetime of the dot states, characterized by the intrin- 
sic line width Tn(w) = 27rgN(w)|VN| 2 and T-p^uj) = 
27T0F,a(w)| Vf\ 2 - For later use, we define Tp(w) = 
[r Fi+ '(u;)+r F _(u>)]/2, which implies P(w)T F (w) = 
[Tp +(<*;) —Tp,- (a;)] /2, or, equivalently, Tp, a (uj) — 
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rp(w)[l + aP(uj)]. Finally, we define the total intrin- 
sic linewidth T(uj) — Tp(uj) + Tn(w). 

The main goal of this paper is to investigate charge 
and spin pumping schemes relying on the variation of the 
ferromagnet's magnetization M(i) = M(t)e p (t), specifi- 
cally its amplitude M(t) or direction e p (t), in time. The 
variation of the magnetization amplitude can be micro- 
scopically modeled by time-dependent shifts of the ma- 
jority and minority bands, E a (t), which leads to a time- 
dependent Stoner splitting AE(t). The variation of the 
magnetization direction, on the other hand, is described 
by time-dependent angles 9{t) and <p(t) in Eq. Q, which 
leads to an explicit time dependence of the operators 

Ofc±(t). 

In order to achieve pumping in the adiabatic regime, 
two system parameters need to be varied in time with 
a relative phase. For a time-dependent magnetization 
direction, these could be the x- and y-component of 
the polarization, P x (t) — P sin 9 (t) cos <p(t) and P y (t) = 
P sin 9 (t) sin ip(t). If, on the other hand, the magnetiza- 
tion direction is fixed and only its amplitude is varied 
via a time-dependent Stoner splitting AE(t), a second 
pumping parameter is needed. Therefore, in the follow- 
ing derivation, we allow for a time-dependent level po- 
sition e (t), which can be experimentally controlled by a 
gate voltage. 



III. FORMALISM 

The problem under consideration is, in general, com- 
plicated since it combines a few interacting degrees of 
freedoms of the quantum dot with a large number of non- 
interacting degrees of freedom in the leads and an explicit 
time dependence. For a time-independent Hamiltonian, 
it is possible to integrate out the lead degrees of freedom 
to obtain an effective description of the system in terms 
of the reduced density matrix. Then one can perform 
a diagrammatic expansion of the time evolution of the 
reduced density ma trix: t o write the kinetic equations of 
the reduced systemP 7 -! 7 -^ 

In the presence of an explicit time dependence of the 
Hamiltonians describing the dot and/or the tunnel cou- 
pling, the kinetic equations can still be formally derived 
but they become much more complicated, such that a full 
solution is only achievable in special cases.^In the limit 
of adiabatic pumping, however, the kinetic equations can 
be simplified considerably by performing an adiabatic ex- 
pansion, i.e., an expansion in the pumping frequency fl 
(or powers in time derivatives of the pumping parame- 
ters), assuming that the response of the system is mu ch 
faster than the change of the system parameters! 42 * 70 ^ 

In the present context, though, also the lead Hamil- 
tonian acquires a time dependence. This possibility is 
not included in the diagrammatic technique presented in 
Ref . l70l where a central step relies on integrating out the 
lead degrees of freedom, for which a time-independent 
equilibrium distribution is assumed. In order to treat 



time-dependent lead Hamiltonians, we start by perform- 
ing an adiabatic expansion already at the very first step 
for the Hamiltonian before integrating out the lead de- 
grees of freedom. The expansion should be performed 
about the time t at which the charge and spin currents 
are calculated. At any different time r, the Hamiltonian 
H{t) = H (t) + V(t) is approximated by 

Ho(t) = i?dot(i) + #lcad N + #lcad f(*) (6) 
V(T) = H tun + (T - t) [fldatW + #lcad F (*)] • (7) 

Higher time derivatives of i?dot and i?i ea d f are neglected 
in the adiabatic expansion. The time derivative of the dot 
Hamiltonian is given by 



H dot (t) = eJ2 4d a . 



(8) 



To derive the time derivative of the ferromagnetic lead, 
we use ctfecr = (in Schrodinger picture) to obtain a, = 

«(^/ 2 ) a L-*(^/2) (a^vu, ka - a ^„u, m j, 
troduced the notation a = —a. This yields 



a cos 9 (V, — sin 9 at w ) , where we in- 



lead F 



(*) 



E<* a la a kc 



+ + Sin ^) a \a a koi , (9) 

which includes non-spin-flip (diagonal) terms in the first 
and spin-flip (off-diagonal) ones on the second line. 

We remark that there is no time derivative of the tun- 
nel Hamiltonian since it does not depend explicitly on 
time, which is evident from Eq. ([5]). For integrating out 
the lead degrees of freedom, however, it is convenient 
to rewrite in the following the tunnelling Hamiltonian 
in terms of the majority /minority electron operators by 
using the time-dependent transformation Eq. Q. 

After having separated the Hamiltonian H (t) into the 
time- independent part Ho(t) of the decoupled system and 
the correction V(t) due to tunneling and the explicit time 
dependence of the system parameters, we can continue 
in a similar way as in the derivation o f the diagrammatic 
rules for a time-independent system.^EU First, we ex- 
press the quantum-statistical expectation value of any 
operator A at time t as integral over the Keldysh con- 
tour K, that runs from time r = — oo to t and then back 
to — oo, 



(A(t)) = tr 



g T K exp 



A" 



drVir)! A(t) 



(10) 



The Keldysh time ordering operator Tk orders all oper- 
ators along the Keldysh contour, the index / indicates 
interaction picture with respect to Ho(t), and po is the 
(full) density matrix at time — oo. The latter is assumed 
to be a tensor product of the density matrices for the 
quantum dot and the leads, with the leads being in an 
equilibrium state determined by H\ ca< \ n and -ffi ca d f(^), 
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respectively. By doing so, we neglect nonequilibrium dis- 
tributions of the ferromagnet. This is justified as long as 
the time scale for spin-relaxation towards equilibrium in 
the ferromagnet is much shorter than the time scale for 
transport, given by H/T. In that case, the non-adiabatic 
corrections of the ferromagnet's density matrix are neg- 
ligible in comparison to the non-adiabatic corrections of 
the transport processes, and we only need to develop a 
theoretical description of the latter. The next step is 
to expand the exponential in powers of V(t)j. In the 
diagrammatic language, each V(t)j is represented by a 
vertex. In the present context, there are three different 
types of vertices: H tun (represented by a full circle) con- 
tains one dot and one lead operator, iTdot(i) (represented 
by an empty circle) two dot operators, and H\ ca aF(t) 
(represented by a double cross) two lead operators. 

At this stage, it is possible to perform the partial trace 
over the lead degrees of freedom. By using Wick's the- 
orem, we contract the lead operators in pairs. In the 
diagrams, each contraction is represented by a tunneling 
line. The diagram for the full time evolution of the re- 
duced density matrix from — oo to t is then a sequence 
of irreducible blocks, defined as the parts of a diagram, 
where a vertical line at any time r crosses at least one 
tunneling line. This infinite series of irreducible blocks 
can be summed up by a Dyson equation. 



itive form 



To obtain the kinetic equation for the matrix elements 
(Xi\Qred(t)\X2) °f the reduced density matrix 
Prod, we use for the operator A in Eq. (10 1, the projector 



|X2><Xi| w it h Xi,X2 e {0,t,l,d}. This leads to 



dt 



P(t) 



dt' W(t,t')p(t') 



(11) 



where p is the vector of all relevant density matrix ele- 
ments, p = (poiPt,Pi.,Pd>p\ ,p[) T - The diagonal matrix 
element p x = p* is the occupation probability of the 

state x °f the dot, while the off-diagonal ones, pj = (pj)* 
describe coherent superpositions of up- and down-spins. 
The kernel W (t, t') is a 6 x 6 matrix which acts on the 
vector p(t'). It represents an irreducible block in the di- 
agrammatic language and describes transitions between 
the states at time t' and time t. 

In order to distinguish the charge and the spin de- 
grees of freedom, it is convenient to parametrize the 
reduced density matrix not by the six elements of the 
vector p but rather use the probability vector P = 
(P ,-Pi,-Prf) T = (po,Pt +PliPd) T and the spin vector 

T ( pj+p\ . p\~p\ p t - P4 , X 



(^1, Sy, S Z ) 



instead. 



jP (t) = J dt' W p (t, t') P (f ) + W s (t, t') S (0 (12a) 

— oo 

i 

S (t) = J dt' M p (t, t') P 0') + M s (t, t') S (t 1 ) . (12b) 



d_ 
dt 



The matrix elements of the 3x3 matrices ~W p (t,tf), 
W s (t,t'), M p (t,t r ), and M s (t,t') are the proper linear 
combinations of the matrix elements of the 6x6 ma- 



trix W (t,f). The first contribution to Eq. (12b I, which 



depends on P, can be interpreted as spin accumulation, 
while the second one describes relaxation and coherent 
rotation of the spin S. 

We remark here that spin rotational symmetry about a 
given axis simplifies the structure of the kinetic equations 
Eqs. (12a I and (12b I. The simplified form of the kinetic 



equations is derived in Appendix [AJ This symmetry ap- 
plies when the ferromagnet's magnetization direction e p 
is constant in time. But even for a time-dependent e p (t), 
the spin rotational symmetry is present for the instanta- 
neous part of the kinetic equations (defined in the next 
section) . 



A. Expansion of the kinetic equation 

Although we have already linearized the time depen- 
dence of the Hamiltonian H(t) about the time t, we still 
need to perform a systematic adiabatic expansion for the 
kinetic equations. For this task, we follow Ref. [70l The 
lowest (instantaneous) order describes the equilibrium 
situation when all parameters are frozen to their values 
at time t, 



o = wf ) p« 



(13) 



The next-order (adiabatic) correction contains all contri- 
butions linear in the pumping frequency i.e., with one 
first-oder time derivative appearing, 



d 
dt 



p W = W (i) p (a) + W W p (i) 



* dt * 



(14) 



Here, the index t indicates the time with respect to which 
the adiabatic expansion has been performed and the su- 
perscript indicates the order of the adiabatic expansion. 

(a) 

By construction, W t contains exactly one empty cir- 
cle or one double cross vertex representing -ffdot(i) or 
-fficadF(i), respectively. Furthermore, we have intro- 
duced the Laplace transform to define the abbreviation 



W 



(i/a) 



w 



(i/a) 



(z) 



/! 0o «ft / W t (i/a: '(t-t / ) as well 



z=0 4 



as <9W t (i) = dW t 



(i) 



(z) /dz 



,0+ 



In this basis, the kinetic equations take the more intu- 



On top of the adiabatic expansion, we perform a sys- 
tematic perturbation expansion in powers of the tunnel- 
coupling strengths T. The order of the expansion in the 



5 



tunnel coupling will be indicated by an integer super- 
script. For example, W^' 1 ^ indicates the first order in T 
of the adiabatic correction of the kernel. This expansion 
is straightforward and can be easily performed (details 
can be found in Ref. 1701 ) 

The perturbation expansion of the kernel starts in first 
order in T, which corresponds to sequential tunneling 
processes described by diagrams containing one tunnel- 
ing line. The perturbation expansion of the instanta- 
neous probability vector starts in zeroth order in T, since 
it corresponds to the time-independent problem with all 
parameters frozen at time t. Thus, the normalization 
conditions read e T P^ 1,0 -* = 1 and e T P^ 1 ' 1 ^ = with 
e T = (1,1, 1) T . The adiabatic correction (i.e. first or- 
der in pumping frequency uS) to the probabilities starts 
in minus first order in T and obeys the normalization con- 
ditions e T p^' -1 ) — and e T p^ a, °' = 0. When going to 
the limit of weak tunnel coupling, T — > 0, one needs to 
keep in mind that the adiabaticity condition, Q <C T, has 
to remain fulfilled. Therefore, p^' -1 ^ oc tt/T does not 
diverge. 



B. Pumped charge and spin current 




FIG. 2. Examples of diagrams belonging to the instanta- 
neous (i) and adiabatic (a) kernel in first order in the tunnel- 
coupling strength F, respectively. 

charge and spin as 

T 

Qx= J lN{t) dt (16a) 
o 

r 

S x = J J N (t)dt, (16b) 
o 

where the index X indicates the pumping parameters. 



We are interested in the pumped charge and the 
pumped spin (projected along a time-independent spin 
quantization axis) through the quantum dot. The 
pumped charge and spin currents flowing into the 
normal-metal lead can be written, respectively, as 

J N (i) = e f dt'e T W N ^ (t, t') p (t') (15a) 
t 

J N (t) = i J dt'e T W N ' S (M')P(O (15b) 

— oo 

where Wf' Q (t, t') = ' T (t, t') + Wf 4 (t, t'), and the 
matrix elements of the kernel 'W t ' a (t,t') are evalu- 
ated with the same diagrammatic rules as for W f (t,t') 
but with the difference that each diagram is multi- 
plied with the number of electrons with spin a enter- 
ing the normal lead during the transition minus the 
ones leaving the normal lead. The z-component of the 
vector W N ' S = (W N < S * , W N < S « , W NA ) is given by 
W n,S; = wf' T (M') - Wf A (t,f), and the x- and y- 
component are obtained by the same expression but with 
the spin quantization axis being chosen along the x- and 
y-direction, respectively. 

In the same way as for the kinetic equation, we perform 
an adiabatic expansion of the charge and spin currents. 
The instantaneous currents vanish since, in the absence 
of a bias voltage the instantaneous terms describe the 
equilibrium situation. 

By integrating the pumped charge and spin currents 



IV. DIAGRAMMATIC RULES 



We explain here how to evaluate the instanta- 



neous diagrams W 



(i) 



Wf (z) 



and <9W 



i) 



z=0 4 



<9W t (i) (z) jdz 



as well as the first adiabatic correc- 



tion W 



(a) 



z=0+ 



over one pumping cycle T 



we obtain the pumped 



A. Rules for the instantaneous diagrams 

In the following, we write down the rules for the instan- 
taneous diagrams w| 1,n '(z) to n-th order in the tunnel 
coupling. An example of a diagram belonging to the in- 
stantaneous kernel element (W^' ^)ot an d its adiabatic 
correction (W| a,1 ')o| are shown in Fig. [2] All diagrams 
contributing to these matrix elements are presented and 
evaluated in Appendix |B| 

(1) Draw all topologically different diagrams with 2n 
full-circle vertices connected in pairs by directed 
tunneling lines. Assign a reservoir index r = N, F, 
an energy u>, and a spin index (<r =t, 4 for r = N 
and a = ± for r = F) to each of these lines. Assign 
quantum-dot states \ £ {0,t,|,d} and the cor- 
responding energies E x e {0, e, e, 2e + U} to each 
element of the Keldysh contour between two ver- 
tices. Furthermore, draw an external line with the 
(imaginary) energy — iz from the upper leftmost be- 
ginning of a dot propagator to the upper rightmost 
end of a dot propagator. 



(> 



(2) For each time segment between two adjacent ver- 
tices (independent on whether they are on the same 
or on opposite branches of the Keldysh contour) 
write a resolvent 1/R, where R is the difference of 
left-going minus right-going energies (including en- 
ergies of tunneling lines and the external line - the 
positive imaginary part of iz will keep all resolvents 
regularized) . 

(3) The contribution of a tunneling line consists of a 
prefactor l/(2w) and a factor Tn(oj) for r = N or 
^F,a{^) for r — F. This is multiplied with f + (uo) — 
f(uj) if the tunneling is going backwards with re- 
spect to the Keldysh contour and f~ (uj) = 1 — f(ui) 
if it is going forward. Here, f(ui) is the Fermi func- 
tion and T F>a (u>) = T F (u>) [1 + a P(u>)] . 

(4) Each full-circle vertex with an incoming or outgo- 
ing tunneling line (with spin a or a for r = N or 
r = F, respectively) changes the dot state assigned 
to the Keldysh contour by adding or removing an 
electron with spin a. In the case r = F, the match- 
ing of the different spin quantization axes of dot and 
lead is achieved introducing a prefactor A aa for a 
vertex with an outgoing line and A aa = (A aa )* for 
a vertex with an incoming line. 

(5) The overall prefactor is given by (—«)(— l) b+c where 
b is the total number of vertices on the backward 
propagator and c the number of crossings of tun- 
neling lines. Furthermore, there is a minus sign for 
each vertex connecting dot states | f) an d |d). 

(6) Integrate over the energies of tunneling lines and 
sum over the reservoirs and the spin indices of the 
leads. 



B. Rules for the adiabatic diagrams 

The adiabatic corrections to the kernels are described 
by diagrams that contain one vertex that is associated 
with the time derivative of the Hamiltonian. This may be 
an empty-circle vertex for (r— t)£/dot(i) or a double-cross 
vertex for (r — £)-Hi ea d f(*)> respectively. No tunneling 
line is attached to the empty-circle vertex since fldot(i) 
does not contain any lead operator. In contrast, there 
are two lead operators in _ffi ca d p(t). As a consequence, 
two tunneling lines are coupled to a double-cross vertex: 
first (with respect to the Keldysh contour) an incoming 
and then an outgoing one. 

When performing Wick's theorem, the two lead oper- 
ators of the double cross may be contracted either with 
each other or with two other lead operators of full-circle 
vertices. The earlier possibility, however, does not con- 
tribute (any diagram with such a self-contracted double 
cross vertex on the upper propagator cancels with the di- 
agram obtained from it by moving the double-cross ver- 
tex to the lower propagator.) 



To evaluate the diagrams W f (z) the rules for 
wf'™' (z) need to be modified in the following way: 

(1') In addition to the 2n full circle vertices there is 
either one open-circle or one double-cross vertex 
which, in principle, can sit everywhere on the con- 
tour. The latter is connected to two tunneling 
lines of the ferromagnet: first (with respect to the 
Keldysh contour) an incoming and then an outgo- 
ing line. The double-cross vertex may be either 
diagonal or off-diagonal in spin. In the first case, 
the two lines carry the same energy uj and the same 
spin a. In the second case, one carries uj and a, and 
the other one uj + aAE and a (= —a). Add an ex- 
ternal frequency line with (imaginary) energy — iz' 
from the empty-circle or double-cross vertex to the 
upper right corner of the diagram. 

(3') When applying rule (3) for the two tunneling lines 
connected to a double-cross vertex, only one prefac- 
tor rF, Q (w)/(27r) (the one of the line which carries 
uj and a) has to be taken into account. 

(4') An empty-circle vertex comes with a factor E x (t), 
where E x {t) is the energy of the dot state x as- 
signed to the Keldysh contour segment where the 
empty circle vertex is placed. For a double-cross 
vertex that is diagonal in spin, the factor is E a (t). 
The factor for a double-cross vertex off-diagonal 

in spin is (AE/2) (^—6 + iaip sin 8^J , where a is the 

spin of the outgoing line. Furthermore, one needs 
to perform a first derivative with respect to z' and 
then send z' to + . 

(5') When applying rule (5), b is the total number of 
all vertices (including full circle, empty circle, and 
double cross vertices). Similarly, the number c in- 
cludes also crossings of tunneling lines at a double 
cross vertex. 

We remark that the contributions of diagrams with the 
rightmost vertex being a double cross cancel out since for 
each such diagram with the double cross on the upper 
propagator, there is a partner diagram with the double 
cross on the lower propagator which only differs by a 
minus sign due to rule (5). 

Furthermore, we remark that for time-independent 
magnetization direction of the ferromagnet, 6 — up — 0, 
only double-cross vertices diagonal in spin space appear. 
In this case, it is possible to directly integrate in time 
over all positions of the empty-circle or double-cross ver- 
tex between two full-circle vertices. As a consequence 
one does not need to explicitly draw this empty-circle or 
double-cross vertex but the adiabatic correction to the 
kernels are obtained from the same diagrams as for the 
instantaneous ones with modified rules. This route has 
been used in Ref. [701 to account for a time dependence of 
the dot level position. 
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V. RESULTS 

We start by deriving the expressions for the charge 
and spin currents to lowest-order in the tunnel coupling 
strength. To this order, only the instantaneous kernels 

(i 1) 

WJ ' to first order in F enter both the instantaneous 
limit and the adiabatic correction of the kinetic equa- 
tions, 



dt 



Pi 





(i,0) 



W (1,D p(.,0) 



W 



(i,l) _(*,-!) 
Pt 



(17) 
(18) 



As a consequence, the kinetic equations and the expres- 
sions for the pumped charge and spin current to this order 
are independent of the pumping scheme, i.e., the choice 
of pumping parameters. Furthermore, there is rotational 
spin symmetry about the ciXIS Cp , 1 . e., the kinetic equa- 
tions take the form derived in Appendix [A] 

After explicit calculations, that are summarized in Ap- 
pendix [Cl we obtain 



4 B,0) (*) 



i-p 2 



r 2 



dt {n 



,(a,0) 



(t) 



P l - 



2T 



l_ P 2 



r 2 



d 

Sp dt { 



,(i,0) 



(19a) 
(19b) 



for the charge and the spin current. The latter is polar- 
ized along e p . Both currents have a similar dependence 
on P, Tn, T f , and the instantaneous part of the aver- 
age number of electrons (n) = Pi + IP^ expanded to 
zcroth order in T. As a consequence, for each moment in 
time, the ratio between the component of the spin cur- 



rent J, 



(a,0) 



X 



(*) 



along the (instantaneous) 



symmetry axis and the charge current is given by 



J, 



(o,0) 



(*) 



,(o,0) 



(*) 



2e T 



(20) 



In order to calculate the pumped charge and spin per 
cycle, we need to specify the pumping scheme. As an- 
nounced in the introduction, we will consider two differ- 
ent scenarios. 



For pumping scheme A, it is convenient to choose the 
same spin quantization axis for the dot and the ferro- 
magnet, + =f and — =],■ The variation of the mag- 
netization amplitude can be microscopically modeled by 
time-dependent Stoner splitting AE = Ek— — Pfc+- 

To be more specific, the majority- and minority-spin 
bands are shifted relative to each other in time, E±(t), 
in such a way that both the Fermi energy and the to- 
tal number of electrons in the ferromagnet remain con- 
stant. This is, at low temperature, fulfilled for E + /E- = 
~ Qf,-(cf)/ Qf,+(cf)- As a consequence, both the polar- 
ization P(t) and tunnel-coupling strength rp(t) vary in 
time. This alone does not establish any adiabatic pump- 
ing since the time variations of P{t) and T F (t) are in 
phase. Therefore, we use the dot-level position e as an- 
other out-of-phase time-dependent parameter. In this 
situation there are two pumping cycles occurring simul- 
taneously: one cycle with {e, P} and another cycle with 
{e,rp} as pumping parameters. 

In the following, we concentrate on the limit of weak 
pumping, i.e. we write the level position e as well as the 
polarization P or the tunnel-coupling strength T F as a 
sum of the average value and a small variation, e(t) = 
l + 5e(t), P(t) = P + SP(t), and T F (t) = f F + ST F (t), 
and expand the pumped charge and spin to bilinear or- 
der in the variations. Then, we integrate the charge 
and spin currents over one pumping cycle to obtain the 
pumped charge and spin, respectively. The latter are 
proportional to the area of the pumping cycle in pa- 
rameter space, quantified by the dimensionless quantities 

Vi = Jo dt 5 P e 5P and V2 = Jo dt 5 fie St F /f F for the 
cycles with {e, P} and {e,r F }, respectively. 

Before addressing the sum of the two pumping con- 
tributions, we analyze them separately. The reason is 
that their relative weight to the total pumped charge and 
spin in the weak-pumping regime depends on the ratio 
v = P 772/771 that, in turn, depends on details of the fer- 
romagnet 's band structure. To be specific it depends, at 
low temperature, on the density of states as well as the 
energy derivative of the density of states of the majority 
and minority spins at the Fermi energy. Expressing them 
in terms of the polarization P = P(e F ) and the total 
density of states g F = g F (e F ) = ^[£>f,+ ( £ f) + £>f,-( £ f)] 



and their derivatives P' = dP(w) / du\ u 
3gF(w)/9w| Ii , =eF leads to 



and 



A. Pumping scheme A: time-dependent 
magnetization amplitude 

In pumping scheme A, we assume that the magnetiza- 
tion amplitude M(t) of the ferromagnet changes in time 
while its directions remains fixed. Experimentally, this 
could be realized by using paramagnetic diluted magnetic 
semiconductors which exhibit a large Zeeman splitting. 
A small, externally applied magnetic field would, then, 
spin polarize the lead, and the degree of spin polarization 
could be varied in time by making the external field time 
dependent. 



p 2 p' 
1-P 2 p 



QF 



_ o P 2 P' 
Z l-P 2 P 



(21) 



For small polarizations P, the ratio v scales as P 2 and is, 
thus, also small. This means that in this case the main 
contribution to pumping is due to the cycle with e and 
P. For large polarizations, P — > 1, on the other hand, 
v — > — 1/2. It is easy to show that for parabolic bands 
the following relation holds — 1/2 < v < 0. 

For arbitrary band structures, however, also positive 
values of v are possible. This can be seen, e.g., for a 
weak ferromagnet with a small Stoner splitting AE by 



expanding the ratio v in AE, 



-AE 2 
4 



0k. 
Qf 



(22) 



which is positive whenever q'^/qf > (q'f/qf) 2 - For ex- 
ample, v > is realized by a functional dependence of 
£>f(w) oc w~ q with a > 0. 



i. Contribution from pumping with e and P 

First, we consider {e,P} as pumping parameters. The 
procedure described above yields for the pumped charge 
Qx and pumped spin Sx — Sx 
axis, 



e p along the symmetry 



P ££ 

± p2 

1 _p2 



r F 

1 T~ 



where r/i = J Q dt 50e SP is the area of the pumping 

cycle in parameter space, (fi)^' ' the instantaneous av- 
erage occupation number, where the level position has 
been replaced by its time average e, and T = Tn + Tp. 

The dependence of the pumped charge and spin on 
the average polarization P and the ratio of the tunnel 
coupling to the ferromagnet and the total coupling, Tp /T, 
is shown in Figs. [3] and |4j respectively. The electrical 
charges and spins flow in different directions, thus the 
particle and spin currents flow in the same direction. 

We find that both the pumped charge and spin vanish 
for T-p/T going to zero, quadratically the former and lin- 
early the latter. For P going to zero the pumped charge 
goes linearly to zero, while the pumped spin remains fi- 
nite. In this limit, the pumping scheme generates a pure 
dc spin current, i.e. a finite spin current with no asso- 
ciated charge current. The spin-pumping efficiency, de- 
fined as the pumped spin per pumped charge, can be 
obtained immediately from Eqs. (23 1 and it reads 



1 + P 



i-P 2 



d 

dpi 
r 2 

2~ 



d 

d/3e 



S (i,0) 



(23a) 



(23b) 



R 



-2e- 



r 



pr T 



(24) 



It is worth noticing that it becomes arbitrarily large for 
Pfp/f -> 0. In the limit Pf F /f ->■ 1 the number of 
pumped charges equals the one of the pumped spins. 



2. Contribution from pumping with e and Fp 

Next, we consider the contribution to the pumped 
charge and spin that originates from pumping with e and 
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FIG. 3. (Color online) Pumped charge in units of 

er\\d$e (n)' 1 ' ' as a function of (a) the time average of the lead 
polarization and (b) the relative tunnel-coupling strength re- 
spectively. The pumping parameters are e and P. 



Tp . This contribution coincides with the results obtained 
for pumping by changing the properties of the scattering 
region (e and Tp) in a F-dot-N structure, investigated in 
Ref. 1421 In the limit of weak pumping, we find 



P 2 rj- rN 
r r 2 r 



z f 



r_F 

r r 



1 -P 2 



pa 



d 

dpi 



v(i,0) 



S f r„ — 



p p 



P 2 n 

1 p2 



O i F 
L f 



l-P 2 



r- 



d 
dpi 



(25a) 



\(i,0) 



(25b) 



where 772 = / dt <5/3e STf /fy is the area of the pumping 
cycle in parameter space, normalized by Tp. The depen- 
dence of the pumped charge and spin as a function of 
P and fp/f is shown in Figs. ^ and As already 
remarked in Ref. 1421 the pumped spin changes signs as a 
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FIG. 4. (Color online) Pumped spin in units of 

— \r\\d$t (n)' 1 ' ^ as a function of (a) the time average of the 
lead polarization and (b) the relative tunnel-coupling strength 
respectively. The pumping parameters are e and P. 



FIG. 5. (Color online) Pumped charge in units of 

e.r\idpe (n)' 1 ' ' as a function of (a) the time average of the lead 
polarization and (b) the relative tunnel-coupling strength re- 
spectively. The pumping parameters are e and Tp. 



function of Tp/P", while the pumped charge does not. 



3. Total pumping efficiency and pure spin current 

After having discussed separately the two contribu- 
tions due to pumping with {e, P} and {e, Tf}, we now 
turn to address the sum of the two. In Fig. ([7]), we show 
the total spin efficiency 



R 



-2c 



5' 



Se,P 



;£,r F 



(26) 



as a function of f f /f for different values of v and a polar- 
ization P = 0.2. For v = (realized for flat bands around 
the Fermi energy) we get only the contribution due to 
pumping with e and P which diverges for T-p/T — > 0. 
A divergence of the spin efficiency indicates a pure spin 
current without a charge current. However, for v = 
the pure spin current is only asymptotically reached for 
Pp/P — > since in this limit the amplitude of the spin 



current vanishes. A negative value of v removes the di- 
vergency (no pure spin current). The most interesting 
case is realized for positive values of v. In this case, the 
divergence of the spin efficiency is shifted to finite values 
of Pf /P, which correspond to a pure spin current of finite 
amplitude. The sign change of R at this point reflects a 
sign change in the total pumped charge current. 



B. Pumping scheme B: rotating lead 
magnet izatiom 

In pumping scheme B, the magnitude M 
of the ferromagnet's magnetization M(i) = 
Me p (t) remains fixed but its direction e p (t) — 

(sin 8 cos ip(t), sin 9 sin </?(i), cos #) T rotates about 
the z-axis. The latter is used as the (time-independent) 
quantization axis for the dot and normal lead electron 
spins a —'f , j., whereas the direction of the majority and 
minority spins a = ± of the ferromagnet changes in 
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time. The Stoner splitting AE — Ek- — Ek+, on the 
other hand, remains constant in time. To experimentally 
induce such a rotation in the ferromagnet, one may 
make use of a ferromagnetic resonance. 



As mentioned above, adiabatic pumping requires the 
time-variation of two system parameters with a rela- 
tive phase. In pumping scheme B, the two parame- 
ters are the x- and y-component of the polarization, 
Px(t) — Ps'md cos(p(t) and P y (t) = P sin 6 sin <p{t). In 
contrast to pumping scheme A, we do not need to vary 
the dot-level position in order to achieve pumping. 



FIG. 6. (Color online) Pumped spin in units of 

— §7?2<9/3£ (n)' 1 ' ^ as a function of (a) the time average of the 
lead polarization and (b) the relative tunnel-coupling strength 
respectively. The pumping parameters are e and IV. 
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FIG. 7. Total efficiency of the spin pump as a function of 
the time averaged ratio Ff /T plotted for different values of v 
when the polarization is chosen to be P — 0.2 . 



In pumping scheme B, the azimuthal angle ip(t) of the 
ferromagnet 's magnetization direction is the only time- 
dependent parameter. As a consequence, the instanta- 
neous avera ge do t occu pation is constant in time. 
From Eqs. (19a) and (19b) we deduce that both the 



charge and the spin current vanish to lowest order in 



the tunnel coupling, J^ a ' (t) 



Oand 3^°\t) = 0. It is, 



therefore, necessary to include the next-order contribu- 
tion in the perturbation expansion in the tunnel-coupling 
strength. 



An explicit calculation, presented in Appendix [D] 



yields a vanishing pumped charge current, L 



(a,l) _' 



0. 



but a finite pumped spin current, which can be nicely 
written in a compact analytical form, 
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)_ r F r N 9 £ (n)^°) 
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' rcl 



rN 
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(Sr/ el ) ; 



x a t e p - 



r F 



B r T s el 



(B rf el ) 



d t e p 



(27) 



Here, and r r ^j are the charge and spin relaxation 
times, respectively, and B describes an interaction- 
induced exchange held that is a consequence of the 
spin-depend ent tun nel coupling of the dot level to the 
ferromagnet! 24 * 25 ^ The explicit expressions are given by 



1 



' rcl 



r[f + (e) + f-(e + U)] 

= r[r(t) + f + (t + u)} 



• rcl 

/>' lw ' //. 



u 



(28) 
(29) 
(30) 



where T = Tn + r F is the total tunnel-coupling strength, 
and f'duj denotes Cauchy's principal value. 

To get the pumped spin per cycle we need to inte- 
grate over one pumping cycle. From the symmetry of 
the problem it is clear that the spin components along 
the x- and y-directions average out, and only the z- 
component survives. Since dte p does not have any z- 
component, it is only the term with e p x <9 t e p that con- 
tributes to the hnite pumped spin per cycle S v . We use 
(e p x dte p )-e z = <p(t) sin 2 and assume a constant angular 
velocity, SI = ip(t), to get 



S v = ]n sin 2 9 G i 1 



(Brf ei y 



Tf 1 + (Br T s el! 



(31) 



expressed with the help of the dimensionless linear con- 
ductance (in units of e 2 /h) 



Gq 



-2tt 



r N r F a e {n)M 



r 2 



(32) 



of the F-dot-N structure for vanishing polarization. 

The exchange held B enters the expression for the 
pumped spin, Eq. (31 1, in combination with the spin re- 



laxation time rf cl . This is quite natural since the latter 
defines the time during which the exchange held acts on 
the quantum-dot spin before it relaxes due to the dot 
electron leaving the dot or another electron entering the 
dot via tunneling to or from the leads, respectively. It 
is remarkable that the presence of the exchange held en- 
hances the pumped spin. To understand this, we notice 
that the exchange held affects the spin dynamics in two 
ways (for a mathematical support of the following phys- 
ical argument, see also Eq. (D6|). It is clear that the 



exchange held should induce a precession of an accumu- 
lated quantum-dot spin about e p (last term in Eq. ( |D6[ )). 
If this were the only effect, the pumped spin into the nor- 
mal lead would read Su- 
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FIG. 8. (Color online) The factor (B r r ^) 2 



/ [l + (B rf cl ) 2 ] 

plotted as a function of the dot level position e for different 
choices of the tunnel couplings, fixed Coulomb energy U = 
20r, and fully polarized ferromagnetic lead (P = 1). 



i.e., the pumped spin would be reducedby this precession. 
There is, however, another contribution in which the ex- 
change held enters, namely the accumulation term (first 
term in Eq. (D6)). This indicates that already during 



the generation of the accumulated spin, the spin dynam- 
ics induced by the exchange held plays a role and gives a 
hnite contribution to the adiabatic correction to the spin 
accumulation. As a result of our calculation, we hnd that 
the combination of the two effects leads to the form in 
Eq. (31 1 for the pumped spin which increases with in- 



creasing exchange held. 

Due to the prefactor r N /r F 



in front of the term 



B 



rel) 



^1 



1+(B 



rcl) 



the exchange held becomes 

more and more important when increasing the tunnel 
coupling to the normal lead. The latter term is plotted 
as a function of the dot level e in Fig. [8] for a constant 
value of the Coulomb on-site energy U = 20T. It goes to 
zero for large values of |e| and has two maxima symmet- 
rically positioned around its minimum at e = -f . With 
increasing Coulomb interaction U, the maxima reach the 
value one and the width of the peaks increases, which 
leads for large values of the Coulomb interaction, as it is 
the case chosen for Fig. [5J to a plateau with a slit. 

In Fig. § we illustrate how the effect of the exchange 
held depend on the tunnel couplings. We begin by choos- 
ing a weak tunnel coupling to the normal lead. 

Then, the second term of the pumped spin including 
the exchange held is small compared to the hrst term and 
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does not play a role. The pumped spin in this situation 
is plotted as a function of the dot level e in panel (a) of 
Fig.[9]for different values of the charging energy U . In the 
absence of Coulomb interaction on the dot, there is one 
peak positioned at e = 0. In the presence of the Coulomb 
interaction on the dot, a second resonance appears at 
e = — U while the amplitude of the maxima is decreased 
and stays constant for Coulomb interactions unequal to 
zero. 

Next, we choose the tunnel couplings symmetrically, 
Tf = Tn- A plot of the spin in this case is shown in 
panel (b) of Fig. [9] The symmetric choice of the tun- 
nel couplings makes the amplitude of the pumped spin 
maximal. One can see that the resonances are still at 
the same positions and that only the amplitude of the 
maxima decreases at first but increases with increasing 
Coulomb interaction. This is due to the enhancement of 
the factor with increasing U. The shape of the factor 
leads as well to a dip at the minimum for U = 10r. 

The pumped spin for asymmetric tunnel couplings 
with a strong coupling to the normal lead is plotted in 
panel (c) of Fig. [9j Compared to the symmetric case, 
the overall amplitude has been reduced and it increases 
more slowly with increasing Coulomb interaction. The 
positions of the main resonances still coincide. For this 
choice of the tunnel couplings the exchange field plays a 
crucial role. The value of the pumped spin is increased 
for gate voltages between the resonances. This gives rise 
to side peaks, one below e = and one above e = — U. 
These peaks overlap with the corresponding main peaks. 
They are visible as individual peaks only for large enough 
values of the Coulomb interaction U. 

We remark that the polarization of the ferromagnet en- 
ters the expression of the pumped spin current (31 ) only 



via the magnitude of the exchange field B. For vanishing 
Stoner splitting AE and, thus, the polarization going to 
zero, Eq. (31) would still give a finite result. However, 



in this case the assumption, discussed in Section III that 
the spin-relaxation time in the ferromagnet needs to be 
longer the the time it takes to complete a transitions be- 
tween minority and majority states is not verified and 
Eq. (31) cannot be applied any longer. When this hap- 



pens the spin flip process in the ferromagnet, described 
by tunneling lines contacted with a double-cross vertex, 
are cut by the spin relaxation time in the ferromagnet. 
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VI. CONCLUSIONS 

We studied adiabatic charge and spin pumping 
through a single-level quantum dot weakly tunnel- 
coupled to a normal and a ferromagnetic lead with time- 
dependent polarization. To this end, we extended a real- 
time diagrammatic approach to account for a time varia- 
tion of the ferromagnet's properties. We investigated two 
different pumping schemes. In the first one, the ampli- 
tude of the ferromagnet's polarization is changed in time. 
To establish pumping, we chose the dot's level position as 



FIG. 9. (Color online) Pumped spin in units of sin 2 9 as a 
function of the level position e of the dot for different values 
of the Coloumb interaction U. The temperature is chosen to 
be ksT — F and the ferromagnetic lead is assumed to be fully 
polarized (P — 1). 

(a) The coupling to the ferromagnetic lead is strong 

(r F = 2or N ). 

(b) The coupling to the leads is chosen symmetric (Tp = IV). 

(c) The coupling to the ferromagnet is weak (20rp = Tn). 
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a second pumping parameter. A pure spin current with- 
out any charge current is only possible for special choices 
of the system parameters. The second pumping scheme 
relies on the rotation of the magnetization direction of 
the ferromagnet. In this case, the pumped charge cur- 
rent always vanishes, i.e., a pure spin current is generated 
without fine tuning of the system parameters. 

We acknowledge financial support from DFG via SPP 
1285 and from EU under Grant No. 238345 (GE- 
OMDISS). 



Appendix A: Kinetic equations in case of rotational 
spin symmetry 

We consider the equations of the from 4? p = 
W p, where the vector p of matrix elements 
of the reduced density matrix is given by p = 
(po,p-f ,pi,pd,Pppl) T ■ Within this appendix, the short- 
cut notation W p could either represent the time convo- 
lution / dt'W(t,t') p(f) or the product W t (i) pf } for 

— oo 

the kinetic equation in the instantaneous limit. 

The basis change from p to P and S is accomplished 
by the transformation 



/l 











\ 


1 


1 











1 
















1/2 


1/2 










-i/2 


i/2 




\0 1/2 


-1/2 








/ 



(Al) 



In case of rotational spin symmetry about the axis e p , it 
is convenient for the following derivation to quantize the 
spin along this symmetry axis. In that basis the kernel 
W reads: 



W 



Woo 


Wot 


w oi 


W od 





\ 


Wfo 




w n 










w lQ 




w u 


w id 








w d0 


w dt 


w dl 


w dd 




















w\\ 





V 














w ll) 



where the zeros for the off diagonal matrix elements in 
the fifth and sixth column and row are a consequence of 
the spin symmetry. 

As a result, the kinetic equations for P and S read 



(S • i p ) , 



d_ 
dt 



S = (v acc • P)i p - 



BSxe 



v - 



where we split the spin vector S = S" +S 1 " into a parallel, 
S" = (S • i p ) i p , and perpendicular part, S 1 - = S — S", 
and we made use of the abbreviations 



/ 



W p = 



W, 



DO 



\ 



t 



W d0 



\Y^w da 



Wo d 

E„w ad 

Wdd 



Wot ~ w oi 



\ 



J 



(w at - w ai ) 
x w dt - w dl 

' \ (W t0 - W l0 ) N 

- = \J2Aw u -w ia ) 

v \ (w u - W id ) j 
\ = - o ( w n - w it + w u - w n ) 



tt 



w. 



•s 



= -^w\\ 



(w\l - w\\) 



SsW 



n 
ii ■ 



The right hand side of kinetic equation for the spin 
is split into three parts. The first one is independent 
of S and describes spin accumulation. The second one 
models the relaxation of the parallel and perpendicular 
components of the accumulated spin. And finally, the 
third term gives rise to a coherent rotation of the spin. 

We observe that the spin symmetry has several conse- 
quences: only the spin component parallel to the sym- 
metry axis i p enters the kinetic equation for P, the spin 
accumulation is along e p , the spin relaxation is rotation- 
ally symmetric about i p , and the spin rotation is about 
the symmetry axis i p . 



Appendix B: Examples of diagrams 



The diagrams contributing to the matrix element 



(W^' 1 ' l )o^ in instantaneous and first order in T are shown 



in Fig. ( 10 ). Applying the diagrammatic rules, we obtain 

/" (w) 



(Wfi)) 



0| 



e — lj + iz 

rN r F (u;) 

2tt 2tt 



=0 + 



[1 + P(u)] cos 2 - 



[l-P( W )]sin 2 - 



/"(e) [r + Pr F cos6»] . 
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-iz -iz 




FIG. 10. Example of the diagrams which are needed to calcu- 
late the matrix element (Wf'^Jof belonging to the instanta- 
neous kernel in first order in the tunnel-coupling strength F. 



We now calculate the adiabatic correction to the in- 
stantaneous kernel in first order in T for finite ip but 8 = 
(this is necessary only for pumping scheme B). Applying 
the rules (l.')-(6.') we have to draw for the kernel el- 
ement (wj n,1 ')oj the diagrams shown in Fig. (JTTJ) and 
get 



( w t Jot = ~y~ sm ' 



E 

*e{+i,-i} 



dui [l + aF(u)] a AE 

Z7T 



d_ 
dz 1 



1 1 

e — cu — a AE + iz e — ui + iz' + iz 



z=0 + 

2 '=o+ 



[/- (w) /- (w + a AE) + /+ («) /- (w + a AS)] 



_9_ 



7? 



^ sin 2 (9 



1 1 

a AE + iz' + iz e — uj + iz' + «z 



T F d e f(e) + — Pr F /"(e) 



^=0+ 
z'=0 + 



[-/- M /+ (u + a AE) + f+ (w) r (w + a AS) 



In the second step, we neglected the energy dependence of P and Tp. It is crucial that before doing this one needs 
to shift the integration variable uj such that AE does not appear anymore explicitly in the integral. For the presented 
example, this is done by 



duj ElM [i + a P( U )] rH 



2tt 



e - uj + i0+ 2ir 



\l + aP] dui 



r H 

e - w + i0+ 



2tt 



e- a; - a AS + i0^ 



, , r F (cj' — a AE) r , , /" (uj') 

du' — — - [1 + aP(w — a AS)] 



/ 



2tt 

duj EiM ^_ aP(w)] r(«) 



e-ui' + i0+ 



2tt 



e - uj + i0+ 

r (w) 



-Ml - a P] / da; 
2vr 1 J / e-a; + £0+ 



Appendix C: Pumped charge and spin current to 
lowest order in tunneling 

We make use of the rotational spin symmetry to 
rewrite the kinetic equations (17) and ( fl8| ), see Ap- 
pendix [3] An explicit calculation of the kernels yields 



= W (i,i) P (i.o) +v (i,i) (Sf' 0) .e p ). 



(Cla) 



o = (v££-p 



,(i,0) 



+ B S (i ^ 0) x e p , (Clb) 



' rcl 



for the instantaneous limit and 



>(i.O) 



dt 



dt 



w (i,i) p(a-i) j. v (i,i) rs( a,_1 ) 



v^ 1 )(Sf- 1, -e p ), (C2a) 



,0,0) _ , , p (a -1) 



)e P 



j(a-l) 



' rcl 



B S^ a " 1} x e p 
(C2b) 



for the adiabatic correction 
The matrix Wp and \ 
kinetic equations for P are given by 



The matrix Wp ' and vector Vp 1 ' 1 ^ appearing in the 
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-iz -iz -iz -iz 




te te 

FIG. 11. (Color online) Diagrams contributing to the matrix element (Wj a,1 ')o| of the adiabatic correction to the kernel in 
first order of the tunnel-coupling strength T. 



( -2/ + (e) 



v^ 1 ) =2Pr 



/ 



rw o 

2/+(e) -/-(e) -f+( e + U) 2f~(e + U) 
f+(e + U) -2f-(e + U) J 

-f-(e) + f + (e + U) 
\ -f + (z + U) J 



(C3) 



(C4) 



In the kinetic equations for the spin we have introduced 



pr, 



(C5) 



as well as the spin relaxation time rf el defined by 1 /rf el 



the interaction-induced 

"/-(") , / + (") 



ex- 



where 



T[f-(e) + f+(e + U)] and 

change field B = ^ =^f^ + 

'pcku denotes Cauchy's principal value. Here and in the 
following, we drop the energy dependence of the tunnel 
coupling T r (oj) = T r and the polarization P(lu) = P. 
The generalization to an arbitrary energy dependence is 
straightforward. 

As discussed in Appendix [A] the right hand side of the 
kinetic equation for the spin, Eq (Clb) has quite an intu- 
itive interpretation. The first term describes spin accu- 
mulation. As a consequence of the spin rotational sym- 



metry, the accumulation is along e p . The second term 
models spin relaxation. For our model, the relaxation 
turns out to be isotropic, i.e., the relaxation times for 
the spin components parallel and perpendicular to the 
symmetry axis e p are identical and denoted by the same 
r rei m t ne following. The third term, describes a coher- 
ent rotation of the accumulated spin about an effective 
exchange field along the magnetization direction of the 
ferromagnetic lead. Its magnitude B depends on the dot 
level position e and is, thus, tunable via the gate volt- 
age. It strongly depends on the Coulomb interaction. 
In fact, when the energy dependence of Tp and P can 
be neglected, the exchange field vanishes in the absence 
of inter action U = 0. In addition to the predicted spin 
rotation) 24 ! 30 ! the exchange field leads to a splitting of 
the Kondo resona ncep^ which has been experimentally 
confirmed recently! 6 ! 17 ! 18 ! 



From the kinetic equations \G\\ and (|C2|), together 

1 and 



with the normalization conditions e T P^ 1 ' ^ 
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e 1 P\ 
spin P 



(a-l) 



we determine the probabilities and the 



(i,0) Tj(a-l) G (i,o) 



and 



(a-l) 



respectively. We 



hnd that the instantaneous probabilities to lowest order 
in the tunnel coupling are the equilibrium values for the 
decoupled dot, thus the spin vanishes, 



j(*,0) _ 



0. 



(C6) 



and the remaining occupation probabilities are simply 
given by Boltzmann factors, 



>(i,o) 



1 



1 + 2e -0« + e-P^+ u ) 



( 



1 

e -f3(2e+U) 



\ 



, (C7) 



We remark that the accumulated spin does not enter the 
expression for the charge current and the probability vec- 
tor does not enter the expression for the spin current. 
This is a consequence of the fact that the currents are 
evaluated in the normal lead, which is not spin polar- 
ized. 

The lowest-order contribution to the adiabatic current, 



given by Eqs. (ClOa) and (ClOb), is linear in f2 and in 



dependent of T, in contrast to the dc current through a 
system with an applied transport voltage, which scales 
with r. Since f2 <C T, the pumped current goes to zero 
for vanishing tunnel coupling as it should. 

Plugging in the results for the probability vector and 
the spin, we finally obtain Eqs. (19a) and (19b I. 



where ft — \/{k^T) is the inverse temperature. It fol- 
lows that the instantaneous average occupation number 



( n )(»-.o) 
order in T is 



_ p(i,0) 



2p( l '°> f the quantum dot in lowest 



(n) m = 2/+(e) 

W f+(e) + f-(e + U) 

The adiabatic corrections are given by 



j(o,-l) 



1 



' rel 



l_ P 2 



p2 



dt * 



(i,0) 



P 1 - 



l_ P 2 



V 2 



(C8) 

(C9a) 
(C9b) 



where is the charge relaxation time given by 1/ = 
T [/ + (e) + /~(e + U)]. The adiabatic correction to the 
probabilities depends on the spin polarization of the fer- 
romagnetic lead and on the tunnel- coupling strengths to 
both leads. Setting the polarization P to zero we obtain 
the result for an N-dot-N structured The accumulated 
spin has only a component along the symmetry axis e p . 

We proceed in the same way for the charge and the 
spin currents flowing into the normal lead. We denote 
the charge current by In and the spin current by Jn- 
The instantaneous currents vanish, since the leads are 
kept at the same chemical potential (no bias voltage is 
applied.) For the adiabatic corrections, we find 



O-i) 



j n W ~ °< r s 

1 r rol 



(ClOa) 
(ClOb) 



where we defined the vector 
/ 



V 



ll-f-(e) + f + (e + U)} 
-f~(e + U) 



Appendix D: Pumped charge and spin current to 
first order in tunneling for pumping scheme B 

In pumping scheme B, the azimuth angle tp(t) is 
the only time-dependent parameter . As a co nsequ ence, 



and S 



(a-i) 



as given by Eqs. ( C9a I and ( C9b I van- 



ish, which leads to a vanishing charge and spin current 
to lowest order in the tunnel coupling, 1^'°^ (t) = and 

J^'°^(£) = 0. It is, therefore, necessary to include the 
next order in the perturbation expansion in the tunnel- 
coupling strength. The adiabatically pumped charge and 
spin currents to next order, 



4 a,1) (*) = 
j^' 1) w = s; 



N( 
(a,0) Tn 



1 

' rel 



(a,0) 



(Dla) 
(Dlb) 



depend on P 



(a,0) 



and 



,o) 



The latter are obtained 



from the kinetic equations to next-to-lowest order in the 
tunnel-coupling strength. 

The kinetic equations expanded to next order 
the tunnel coupling simplify if we use P 



(a,- 







and S 



(i,0) 



s; 



(a-l) 



0, resulting from Eqs. (C9a 



and ( C9b ) . Furthermore, when solving the higher-order 



kinetic equations, it turns out that Pt'' 1 ^ is constant in 



time and that P 



(a,o) _ 



0. This can be easily understood 



noting that the rotation of the fcrromagnet's magnetiza- 
tion direction does not affect the probability distribution 
for empty, single, and double occupation. It immediately 
follows that the pumped charge current is always zero. 
Furthermore, we find that ^ is always parallel to e p . 
This is consistent with the fact that in the instantaneous 
limit e p is a symmetry axis. 

To simplify the presentation, we immediately make use 
of these results when writing down and solving the kinetic 
equations in the following. To obtain the pumped spin 
current, we need S^ - 1 . The latter is determined from 
the adiabatic correction to the kinetic equation for the 
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spin, 



and multiply it from the left with q T = (0,1,2), since 



W 



(U) 



2r/(pr F )(v 



(i,l)\T 



and q 3 



dt 



S (U) = M (a,l) p(i,0) 



;(a,0) 



,(1,1) 
j(i,l) 



B S 



(a,0) 



x e v 



(D2) 



-2Pr F /(rr r 6 el ). This allows us to solve for Sf ' and 
arrive at 



j(i,l) 



' rcl 



(i,2) p (i,0) 
v acc * ^ f 



2T 



^Wp'^Pj' 



All other terms that would formally appear in the expan- 
sion are vanishing, as mentioned above. While rf el and B 

(a 1) 

are already known, and M p ' is straightforwardly con- 
structed by applying the diagrammatic rules, the spin the diagrammatic rules 



l-P 2 



(D5) 



V 2 



Now, vice'' and w{, 1,2 ' can be evaluated with the help of 



s[ entering on the left hand side is still unknown. To 
determine the latter, we write down the instantaneous 
kinetic equations for the spin in next-to-lowest order in 
the tunnel coupling, which immediately yields 



s; 



(i,i) 



' rcl 



\ acc 



p 



0,1) 



,(1,2) 



•Pi 



(1,0) 



(D3) 



We remark that the kinetic equation ( D2 1 for the adi- 



abatic correction s[ a '°' > to the spin contains two source 
terms. One is given by the time derivative of the in- 
stantaneous spin Sj . Since the latter is, for symmetry 
reasons, directed along e p , the time derivative is along 
dti v . The other source term is the adiabatic correction 



to the spin accumulation Mp 
nents along e p x <9 t e p and dti p 



(a,l) p (i,0) 



which has compo- 



Collccting everything, Eq. (D2| reads 



but depends on Pj . To close the set of equations, we 
write down the instantaneous kinetic equations for the 
probabilities in next-to-lowest order in the tunnel cou- 
pling, 



= k lye p x d t e p 



B T? el d t e p 



j(a,0) 



' rcl 



B S 



(a,0) 



x e E 



(D6) 



= w^ 1 ) Pf' 1} 



W (i,2) p(l,0) 

* T P t 



(i,i) 



(D4) 



and we can eventually solve for the adiabatic correction of 



the spin and plug this into Eq. (Dlbl to get the pumped 
spin current as given in Eq. (27). 
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